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PREFACE 

This report describes part of a comprehensive, continuing research program 
in remote sensing of the environment by aircraft and satellite. The research is 
being conducted for NASA's Lyndon F Johnson Space Center, Houston, Texas, by 
the Environmental Research Institute of Michigan (ERIM), formerly the Willow 
Run Laboratories, Institute of Science and Technology, The University of Michigan. 
The main objective of this program is to develop remote sensing as a practical 
tool for obtaining extensive environmental information quickly and economically. 

Remote sensing of the environment involves the transfer of radiation from 
objects of interest on the Earth's surface through the atmosphere to a sensor 
either within or outside the atmosphere. Since this intervening atmospheric me- 
dium alters the intrinsic radiation from surface elements, the interpretation of 
multispectral scanner data can be significantly affected. This report deals with 
the definition of the natural radiation field encountered in a realistic remote sens- 
ing environment from the viewpoints of radiative-transfer theory and atmospheric 
optics. 

The research described in this report vas performed under Contract NAS9- 
9784, Task IV, and covers the period from 1 February 1973 through 31 October 
1973. Dr. Andrew Potter served as Technical Monitor. The program was 
directed by R. R. Legault, Vice President of ERIM, J. D, Erickson, Principal In- 
vestigacor and Head of ERIM's Information Systems and Analysis Department, 
and R. F. Nalepka, Head of ERIM's Multispectral Analysis Section. The ERIM 
number of this report is 190100-24-T. 

The author wishes to acknowledge the direction provided by R. R. Legault and 
J. D. Erickson. Helpful suggestions were made by R. F. Nalepka. Computer pro- 
gramming and other technical analysis were provided by N. A. Contaxes, W. B. 
Morgan, and L. R. Ziegler. The author also thanks R. M. Coleman for her secre- 
tarial assistance in the preparation of this report. 
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SUMMARY 

As part of the continuing development of an atmospheric -radiative -transfer model which 
can be used to correct multispectral data, we must consider as many atmospheric effects as 
possible that may have ?n influence on the data. In previous studies we have constructed a 
radiative-transfer model which Included the effect of multiple scattering in homogeneous, non- 
absorbing, plane -parallel atmospheres bounded by spatially uniform surface conditions. 

Although variable atmospheric states characterized by visual range were included, we consider- 
ed only non -absorbing aerosol particles in the aw*nOsphere, Also, the surface dealt with was a 
spatially uniform Lambertian one. 

In the present further work we consider a variety of more realistic atmospheric conditions 
characterized by visual range, composition of aerosols, size distribution of aerosol pi^rticles, 
and ozone absorption. Hence, we have expanded our atmospheric -radiative-transfer model to 
include almost every kind of atmospheric state possible wiMi regard to the composition and 
structure of aerosols. Excluded are certain strong gaseou ; absorption effects in the infrared, 
with this omission justified on the grounds that in the design oensors, one usually ignores 
these regions anyway. 

Another problem addressed in the current study the occasional influence of background 
on the target through radiation being scattered from element .-v outside the instantaneous field 
of view into that field of view. In this report we examine this effect and employ the necessary 
approximations to solve the problem analytically Also, a solution is presented for the case of 
an airborne or spacecraft sensor viewing a non- Lambertian surface from any point in or out- 
side of the atmosphere. This solution can be used in an analysis of the changes to be expected 
in multispectral data as a sensor moves in a horizontal plane over the terrain being investigated. 
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INTRODUCTION 

In the remote sensing of terrain features on Earth^s surf a. e, discrimination of specific 
targets from their backgrounds is essential. One method of accomplishing this is to analyze 
the radiation received in several wavelength intervals on the assumption that classes of objects 
can be defined In terms of their spectral properties. Such analysis depends upon variabilities 
in remotely sensed data which can be attributed to the complex spectral properties intrinsic to 
materials found on the Earth's surface. Discrimination techniques are by no means limited to 
spectral variations alone; other schemes car. be based upon polarization or goniometric character- 
istics of various classes of objects to provide a "signature" and thus a unique means of identifica- 
tion. But remote sensing based upon polarization and/or goniometric properties has not found 
widespread use, probably because (1) very little is known about the intrinsic polarization char- 
acteristics of natural materials and even less concerning their goniometric properties, and (2) 
the sensor requirements and operational systems needed to acquire meaningful data based upon 
these properties are rather complex. In this report, however, we deal exclusively w*th aiiai tical 
approaches applicable to the development and implementation of discrimination tech.aques based 
upon spectral variations. 

Aside from the intrinsic properties of natural materials, there exist variations, such as the 
variability of the atmosphere, resulting from elements extrinsic to these materials. Although 
one would like to collect data under ideal environmental conditions, the implementation of an 
operational remote sensing system necessitates the acquisition of data under conditions which, 
unfortunately, are oftentimes far from ideal. For this reason it is important that the more sys- 
tematic variations associated with atmospheric effects be considered in the analysis of remote 
sensor data. 

The spectral radiance at a sensor is given by the following: 

Lt = LgT . Lp 

where Lg is the spectral radlaiice at the surface 

T is the spectral transmittance between the sensor and the object being observed 
is the spectral path radiance — that is, the radiation along the path between the sensor 
and target arising from multiple scattering or intrinsic emission 

All three of these quantities depend upon the state of the atmosphere; they can vary considerably 
with altitude, wavelength, view angle, sun angle, and surface albedo. A complete atmospheric 
radiation model should provide those quantities which determine a unique set of parameters 
appropriate to the environmental conditions which exist at a particular time. 
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Previous work on this problem has resulted in an atmospheric -radiative -transfer model 
which accov'nts — assuming multiple scattering in cloudless, plane -parallel atmospheres with ^ 

no absorption — for variations in the radiometric quantities named above. Tli i usual assump- 
tion in atmospheric -radiative -transfer problems is that there is little absorption by aerosols 
in the visible part of the electromagnetic spectrum. This is definitely not the case, hcwever, ^ 

■j 

for the contaminated air commonly found in large urban areas where sulfate, carbonate, and 

other soot particles exist in large quantities. Furthermore, from recent investigations of the : 

atmospheric particul ^te load around the world, it is becoming evident that aerosols play an ; 

important role in determining the overall energy balance of Earth. For this reason we inves- 
tigated the influence of contaminated atmospheres on the natural radiation field. 

Another effect treated in this report is the influence of background on target via scattering 
in the atmosphere. Radiation reflected from background objects will eventually be scattered 
into the field of view of a sensor and, hence, will influence the signal from the target. This 
complicated effect requires that the radiative -transfer equation be solved in two or three 
dimensions in order to define the radiation field anywhere in the atmosphere for non-bomoge- 
neous boundaries. In our study of this phenomenon, we assume it to be a function of basic pa- 
rameters describing pertinent «)ptical characteristics of the atmosphere. 
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THE ATMOSPHERIC STATE 

In the remote sensing of Earth's surface by means of electromagnetic radiation reflected 
from terrain objects, the atmosphere can scatter and absorb f' r?.diation and, hence, alter 
the intrinsic radiation characteristics of the object being investigated. The resultant rad.^tion 
received by a sensor depends upon factors such as viewing geometry, time of year, surface 
conditions, and atmospheric state. It is the state of the atmosphere that we shall discuss in this 
section. 

3.1 COMPOSITION OF THE ATMOSPHERE 

The state of the atmosphere at any location and at any time can be defined by a unique set 
of parameters fully describing all aspects of the atmosphere. However, we are not interested 
in all aspects of the atmosphere, but only in those of key importance to remote sensing applica- 
tions. Thus, micro-fluctuations in density or ion content are of no import, but we are certainly 
interested in visibility. Our goal, therefore, is a somewhat limited specification of parameters 
which can be used to define Ihe state of the atmosphere. For purposes of the present study, 
we shall be concerned mainly with the scattering and absorbing properties of gases and partic 
ulates that exist in our atmosphere. 

3.1.1 ATMOSPHERIC GASES 

By definition, the visible portion of the electromagnetic spectrum is that part of it to vhi^h 
the human eye is most sensitive; thus it is not surprising that this part of the spectrum under- 
goes little absorption by the atmosphere. The major permanent gaseous components of our 
atmosphere are molecular oxygen, nitrogen, and argon. Th' ^c gases absorb very little radia- 
tion in the visible and near -inf rared part of the spectrum; only variable components such as 
water vapor, corbon dioxide, ozone, sulfur dioxide, nitrogen compounds, methane, and other 
trace gases absorb juch radiation to any extent. Of these variable components, ozone, water 
vapor, and carbon dioxide are the most important absorbers of radiation in the near -infrared. 
For the purposes of v‘his study, we shall consider only ozone in the Chappuls bands in the 
region 0.44 to 0,74 /iir* and assume that all other radiation in the region 0.40 to *-2.0 fxm is not 
absorbed by gases. Since the gaseous absorption regimes of the various atmospheric compo- 
nents are well knowTi, it is relatively easy to select those spectral regions within which absorp- 
tion is at a minimum. More detailed studies of atmospheric gases have been done by Goody 
[1| and Zuyev (2|. 

Throughout the visible part of the spectrum, ozone is the primary gaseous absorber. 

There are variable amounts in the troposphere, but the general altitude variation is known [3 1 . 

A profile of the ozone number density is ih tsir ^ed in Fig, 1. Note that the maximum density 
occurs at an altitude of about 23 km. 
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3.1.2 AEROSOLS 

One mc.y define an aerosol as a semi -permanent suspension of solid or liquid particles in 
the Earth’s atmosphere. Typical aerosol distributions would be hazes, clouds, mists, fogs, 
smc^es, and dusts. It is the aerosol component of the atmosphere which is quite variable. 
Depending upon location and time, aerosols can have many compositions, sizes, shapes, and 
densities. In condensation processes, the predominant composition is water (for which the re- 
fractive index is well known). Depending upon the history of the distribution, however, various 
contaminants can mix with water and water vapor in the formation of droplets to produce a 
composition quite different from that of water. The major primary source of atmospheric 
aerosols on a world-wide basis is sea spray [4]; other primary sources are wind-blown dust 
and smoke from forest fires. While the compositioi of worldwide aerosol distribution can be 
estimated, this composition can nevertheless vary over a wide range depending upon its prox- 
imity to various local sources. The only ways m which the composition can be determined is 
by direct m situ sampling of the air in a particular location or by radiometric techniques. 
Much in situ sampling has been done by Volz [5, 6], Grams et al. [7|, and Flanigan and 
Delong [8]. As a result of their work and that of others, the mean world-wide aerosol can be 
considered to have a real refractive index of about l.b and an imaginary part between 0.01 and 
1.0. Kence, 

m(A) = m^ - img (1) 

where m ^ and m 2 , respectively, are the real and imaginary part of the refractive index, and 
A is the wavelength. For realistic atmospheric conditions, 

1.33 <m^< 1.55 (?) 

and 

0 < m 2 ? 1.0 (3) 

which are roughly independent of wavelength. 

particles also come in various sizes. Generally, however, they can be put into three 

-7 -5 

categories: the Aitken nuclei with radii between 10 and lO cm, the large particles with 

-5 -4 -4 

radii between 10 and 10 cm, and the so-called giant particles with radii greater than 10 

cm. In most hazes, the optically active region made up of those particles in the large or 

giant categories. Junge [9| showed that most aerosol distributions follow the simple power 

law. 

N(z, r) - Cfz)r ^ (4) 
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where N(z, r) is the particle number density for radius, r, at altitude, z 
If is the exponent of tne power law 

Experimentally, v has been found to vary f . i 2 to 5 for various tropospheric aerosol distribu- 
tions. Also, other investigators have found a better fit to the particulate d^ta by using the mod- 
ified gamma distrioution, 

N(z, r) = ar^ exp (-br^); 0 5 r < (5) 

where a, a, b, and y are parameters which describe the distribution. 

Deirmendjian [lOj has used this function to characterize different hazes. The following is 
a list of the hazes Deirmendjian considered and the relevant parameters. 


Haze Type a 

M 5.3333 X 10“^ 

L 4.9757 X 10® 

H 4.0000 X 10^ 


a 

r 

b 

1 

1/2 

8.9443 

2 

1/2 

15.1186 

2 

1 

20.0000 


Haze M, used to describe a marine or coastal haze, has a peak in the distribution at 
r = 0.05/um. The haze L represents a continental distribution and has a peak at r = 0.07^m. 

The haze H model can be used to represent a stratospheric aerosol or dust layers: it has a 
peak at 0.10/im A graph of the three hazes is illustrated in Fig. 2. We shall make use of these 
hazes in the detailed analysis of atmospheric radiation. 

For the liquid aerosols it can be assumed that the particles are spherical or nearly spher- 
ical in shape. For solid particles, however, the shape may assume any form. A number of 
investigators [11, 12, 13] have studied the influence of particle shape on the scattering of radia- 
tion. The scattering of electromagnetic l adiation by odd-shaped particles has a different pat- 
tern from that produced by spherical particles; but given a polydlsperse collection of odd- 
shaped particles, the nature of the difference is not clear. Most of the work on radiation in 
atmospheres concerns spherical aerosols; we shall follow suit and neglect the complications 
of panicle jhape. 


Thus far :e have considered the composition, sizes, and shapes of aerosol particles, but 
in order to define the atmospheric state we must also know the concentration of particles. 
Wiegand [14] was the first to measure the vertical profile of an aerosol number concentration 
of condensation nuclei. As a result of his measurements and many others cited in Ivlev [15 1 , 
it was found that the concentration of condensation nuclei obeys an exponential law with altitude. 
It was also determined that a zone of increased concentration of large particles exists in the 
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FIGURE 2. HAZE -TYPE DISTRIBUTION 
FUNCTIONS USED. Units depend on the 
particular model. [10] 
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17 to 23 km altitude range (Junge layer) and in a probable layer under the tropopause at 9 to 10 
km. rhese layers are relatively stable compared to the lower part of the troposphere. As an 
example of an aerosol model atmosphere, Zuyev [2] has constructed the following alt’-' ide- 
dependent, number density-dis^ribution-function for aerosols: 
fN(0)e**^^; 2 ^ 5 km 


N(z) = 


0.03: 5 km ^ z < 15 km 
0.03e^‘^®^, 15 km < z ^ 20 km 
O.Ole'® ®®*: z 5 20 km 


where N(z) and N(0), respectively, are the number densities at altitude z and the Earth’s surface 
(altitude 0). 

In conclusion, we can say that the most significant aspect of aerosol particles in the atmo- 
sphere is their high degree of variability — in composition, size distribution, and especially in 
number density. All the modeLs in the current literature deal with a highly approximate 
average-condition from which large deviations can occur in real situations. The mam^ details 
of aerosol science will not be considered in this report. For a more complete study of the 
physics and chemistry of aerosols, see Mason [16| , Fuchs [17] , Davies [18], or Green and 
Lane [19). 


3.2 TURBIDITY 

Having examined the basic characteristics of the gases and aerosol particles composing 
the atmosphere, we now consider those parameters which relate to the attenuation of radiation 
passing through the atmosphere. 

If a collimated beam of monochromatic radiation is Incident upon a scattering and absorb- 
ing medium, then the intensity of the beam at distance x is given by 

I(x) = (7) 

where I(o) is the beam intensity at the origin. The quantity k , called the volume extinction 
coefficient, is equal to the sum of the volume absorption and volume scattering coefficients. 
Since attenuation of radiation in the atmosphere is caused primarily by molecular scattering, 
scattering and absorption by water droplets, and also by dust particles, Llnke [20] proposed 
that the attenuation be measured by taking the ratio of the total Integrated attenuation coeffi- 
cients to the pure Rayleigh integrated coefficients. Thus, 
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t(A) = 


CO rr OO 

J/Cj^(A)dz + jK^(A)dz + j«jj(A)dz 


OP 


.(x)dz 


( 8 ) 


where k j^(A), and are the volume extinction coefficients for Rayleigh, water aerosol, 

and dust particle scattering respectively. Thus, the turbidity, t(A), is a measure of the depar- 
ture of a real atmosphere from tiie ideal pure Rayleigh atmosphere. Equation (8) can also be 
written as 


t = 1 + W + R (9) 

where W is the humid turbidity factor and R is the residual turbidity factor, Fo* a pure 
spliere free of water and dust, t = 1. Values of the turbidity vary from 3,59 for a continental 
tropical air mass in the summer months to 2.16 for a sea arctic air mass in winter. Tables 
of typical turbidity values are given by Kondratyev [21] for many locations and weather condi- 
tions. The turbidity as a function of altitude is shown in Fig. 3 [22 j ; this figure was constructed 
from data collected by Elterman [23] in optical searchlight measurements. The maximum near 
the tropopause is the result of convective activity causing particles to concentrate at that stable 
position. 


3.3 VISIBILITY 

Visibility in meteorology refers to the transparency of the atmosphere to visible radiation. 
Usually characterized by a quantity, V, called visual range, it depends upon the optical prop- 
erties of the atmosphere, the properties of both the object sighted and its background, and 
illumination conditions. Middleton [24] derives the so-called air-light equation given by 

L = L^e + Lp(l-e ) (10) 

where L is the intrinsic radiance of the ct 
o 

L is the radiance along an infinite path through the atmosphere 
K is the volume extinction coefficient 


Tverskol [25] derives a contrast equation: 

c = 


O 


( 11 ) 


where Is the ratio of the difference between background and intrinsic radiances to background 
radiances, and b^is a brightness factor. This equation can be simplified: 
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C = (12) 

The visual range in km is that distance at which the relative contrast is 2%, i.e., 

f n(C/Cg) = C n 0.02 = (13) 

or 

„ (n 0.02 3.912 

v = = (14) 

Equation (14) is usually taken as the defining equation for visual range at a wavelength of 
0.55/ini, which is near the peak of the human visual-response curve. 
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4 

OPTICAL PARAMETERS 

Section D considered the physical state of the atmosphere. In this section we deal with 
the optical properties of ihe atmosphere; specifically, we define and calculate the appropriate 
atmospheric optical parameters used in radiative-transfer models. 

4.1 SCATTERING THEORY 

Radiation which passes through a homogeneous medium free of any discontinuities will 
cause waves to interfere in such a way that no scattering takes place. Realistic medir, at tem- 
peratures above absolute zero, however, have discontinuities such as crystal impurities, den- 
sity fluctuations, and various sized particles; Eartn’s atmosphere is such a medium. We 
assume, moreover, that all scatterings are independent, that is, that there is no phase relation 
between scattered waves. For this to be true, the distance between scattering centers must 
be at least several times the size of the scattering centers themselves. This condition is 
certainly fulfilled for atmospheric hazes and even for dense fogs. Given this r'“vsical setting, 
we will now determine the basic optical quantities needed for radiative -transfer analysis. 


4.1.1 RAYLEIGH SCATTERING 

The theory of light scattering in the atmosphere originated when Lord Rayleigh [26] found 
an explanation for the bine color of the clear sky. It is nov known that this hue is the result of 
the scattering of radiation from density fluctuations, rather than from molecules as had been 
assumed earlier. If we hive a plane electromagnetic wave with electric field strength 

(15) 

where k is the complex propagation vector 
r is a coordinate position vector 
0 ) is the angular frequency 

then this wave impinging upon a scattering centei whose size is much smaller than the wave- 
length cf the radiation will produce scattered waves in an effect called Rayleigh scattering. 

The scattering cross-section can be obtained by integrating the Poynting vector over a period 
of oscillation* The result is 




Bw\^- 


if 


6 + 36 


( 16 ) 
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where N is the molecular number density 
\ is the wavelength 

6 is an anisotropy parameter to account for polarization effects 
m is the refractive index of the medium 

Of importance here is the 1/X dependence of the cross-section, a characteristic of radiation 
from an oscillating dipole field. 

It can also be shown that the scattering pattern is described by the following formula: 
p(ccs X) = |(1 + cos^ X) (17) 


in which x is the angle between the incoming direction and the scattered direction. The funcfion 
p(cos x), which is called the scattering phase function, is normalized to unity over all 4 t : ster- 
adians, that is, 


jp(cos x)dQ = 1 
Q 


(?8) 


where SI is the solid angle. 


4.1.2 MIE SCATTERING 


The scattering of electromagnetic radiation by homogeneous dielectric spheres is called 
Mie scattering. In principle, determining the scattering and absorptiqp cross-sections anc* 
the scattering phase function is straightforward. First, the wave equation is solved inside 
and outside the sphere; then the solutions are matched at the boundary to determine the con- 
stants. The intensity of the radiation scattered into an angle x is given by 


I - 


^ 2^2 
47t R 



(19) 


where R is the distance from the sphere, and I is the incident intensity. The scattering 

2 ^ 
amplitudes js^ r and js., j' are given by 


_ 2f + 1 r . , 

«=1 

„ ^ 2f + 1 r u _ 1 

^^2 "" V J 

?=i 

In Eq. (20) a^ and b^, which are the Mie coefficients, are given by 
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2pl 


^ ^ (mx)tf/ ^(x) - mtl/ ^ (x) 

^ ■ 4/^(mx)?^(x) - nu//jj(mx)Cj(x) 


b = 


nu//’(mx)i//j(x) - \l/^(mx)\l/'i^(x) 
m\}/'^(mxK ^(x) - <;/j(mxK^(x) 


( 21 ) 


( 22 ) 


where m and x, respectively, are the complex refractive index and size parameters — that is, 
m = m^ - im^ and x = 27rr/X when r is the particle radius. The functions and are the 
Riccati-Bessel functions, and a prime indicates differentiatron with respect to the argument. 
The and functions are given by 

dP-(cos x) 




(23) 


and 


dir .(cos x) 

Tj(cos <) = cos X)Tj(cos x) - sin 


where Pj(co3 x) is a Legendre polynomial. The scattering cross-section is 

. 2 00 


2 ^ 
r, m) = irr Q^(x, m) = g- 


1) 


2 . 2 ' 


(24) 


(25) 


where ©^(x, m) is called the scattering efficiency factor. The total (scattering plus absorption) 
cross-section is 


, 2 00 


a^(\. r, m) = 7rr^Q^(x, m) =^y^(2C + l)Re(a^ + b^ 


(26) 


where Q^(x, m) is the total efficiency factor. Likewise, the absorption cross-section is then 
a^(X. r. m) = a^(X, r, m) - a^(X, r, m) (27) 

with a corresponding efficiency factor Q^(x, m). As in the Rayleigh case, a scattering phase 

a 

function can be defined: 

1 


P(cos x) = 


2irx Qg(x, m) 




(28) 


Calculating all the tunctions given above for various refractive indices and size param- 
eters is quite involved. Nevertheless, computer programs have been written which allow the 
performance o this analysis. Having obtained a program written by J. V. Dave, we used it to 
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calculate the cross-sections and phase functions. For example, we have calculated the scatter- 
ing efficiency factor for homogeneous spheres of refractive i^'dices m = 1.29, 1.29 - 0.04651, 
and 1.28 - l,37i. The results are shown in Fig. 4, It should be noted that the efficiency factor 
is greatest for the real index (that is, when there is no absorption) and decreases with increasing 
imaginary index. The absorption efficiency factor has also been calculated for m - 1.28 - 0.0465i 
and 1.28 - 1.37i. Here the efficiency factor is large for a high imaginary index and small parti- 
cles, as illustrated in Fig. 5. Finally, the total efficiency factor is shown in Fig. 6 for the same 
set of refractive indices. Thus, it caa be seen that efficiency factors vary strongly for different 
refractive indices and size parameters. 


4.2 ATTENUATION COEFFICIENTS 

Knowing the cross-sections, one can then calculate the scattering, absorption, and extinc- 
tion coefficients by multiplying the cross-sections by the particle number density. For real 
atmospheric conditions characterized by haze, togs, and dusts, there is a distribution of particle 
sizes. (One distribution, characterized in Section 3,1.2, is the modified gamma distribution.) 
Thus, for a polydispersion, one must integrate over particle size to obtain the absorption, 
scattering, and extinction coefficients: 

00 

a^(\. m. z) = r)n(z, r)dr (29) 


00 

m, z) = m, r)n(z, r)dr 

0 ■’ 


(30) 


X 

K^(x. m. z) = f m. r)n(z, r)dr 

0 


(31) 


where a, and k denote absorption, scattering, and extinction coefficients 
A designates aerosol 

n(z, r) is the aerosol -particle number density at altitude z for particles in size range 
Ar at size r 

The number density is normalized as follows: 


N(z) = Jii(z, r)dr (32) 

0 

where n(z) Is just the total p**;uoer density. Lllcewlse, the corresponding coefficients for mo- 
lecular scattering o'' .i oe found, The scattering cross-section la given by Eq. (16), and the ab- 
sorption rr^ris-section is usually taken to be that for ozone in the visible spectral region. Thus, 
we have for the complete atmosphere: 






SCATTERING EFFICIENCY FACTOR (Q„) 
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SIZE PARAMETER (x = 2Tr/\) 

FIGURE 4. SCATTEFONG EFFICIENCY FACrOK FOR HOMOGENEOUS 
SPHERES OF COMPLEX REFRACTIVE INDEX m 
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FIGURE 5. ABSORPTION EFFICIENCY FACTOR FOR HOMOGENEOUS SPHERES 
OF COMPLEX REFRACTIVE INDEX m 
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a(X, m. s, z) = z ) + ®A^^* 

(33) 

pfX, m, s, z) = z) ni, s, z) 

(34) 

x(X, s, z) = z) + 

(35) 


where we have shown the explicit dependence on a complex refractive index m and ? particular 
size distribution s. 

It is sometimes useful to deal with an average cross-section for aerosols. Assuming 
that the size distribution is altitude -independent, we can write Eqs. (29), (30), and (31) as 


m. z) = N(z)a^^^(\, m) (36) 

m. 2) = g m) (37) 

K^(X, m, z) = N(z)a^^^(\, tn) (38) 

where 

00 

a ,^^(X, m) = jff. ^(X, m, r)\l/(r)dT (39) 


in which the function is normalized to 1 auu the index i can indicate a, s, or t. We have 
computed many values of the average cross-sections for homogeneous spheres of various re- 
fractive indices and size distributions. The results of this analysis will be presented later. 

4.3 OPTICAL DEPTH 

For radiative-transfer calculations it is usually more meaningful to deal with the dimen- 
sionless quantity, optical depth in a medium, rather than with actual distances. We define 
optical depth as 

t(X, h) = I k(X, z)dz (40) 

h 

and optical thickness as 

X 

Tq(X) = |k(X, z)dz (41) 

where h is some definite altitude. Thus, at the top of the atmos?phere r(A , h-^x) = 0. while at 
the bottom t(\, h=0) = t^(x). Optical depth can be thought of as the distance into a medium 
expressed in units of mean free photon paths. A small Indlcatet^ that little attenuation takes 
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place, whereas a large means that the atmosphere is either strongly absorbing or scatters 
much of the radiation. Optical depths can easily be obtained for Rayloigh atmospheres. 
Elterman [27] has tabulated the results for all altitudes from 0 to 50 km and for selected wave- 
lengths from 0.27 ijxn to 4.00 pm. The optical thickness varies from 1.928 at 0.27 pm to 0.001 
at 1.67 pm. 

Likewise, if the extinction coefficient for aerosols is available, the corresponding aerosol 
optical depths can be determined. By analyzing many experimentally determined aerosol pro- 
files, Elterman [28] has calculated the extinction coefficients and optical depths for realistic 
atmospheric conditions. For a 2-km visual range, the aerosol optical thickness is 2.521 at 
0.36 pm and 1.053 at 0.90 pm. These results indicate that multiple scauering occurs in Earth's 
atmosphere since the mean free photon paths are short compared to the actual distances 
traveled. 


4.4 SINGLE-SCATTERING ALBEDO 

A very important parameter in radiative-transfer analysis is the single -scattering albedo, 
defined as 


^p(X, z) + (3 . in, s, z) 
m, s, z) = K(\, m, s, z) 


(42) 


where k (X, m, s, z) is the total (Rayleigh pli jrosol) extinction coefficient. The albedo 
o) (X, m, s, z) is the fraction of scattering which can occur. Thus, if there is neither aerosol 
absorption nor ozone absorption, then w^(X, m, s, z) = 1 and we have a pure scattering atmo- 
sphere. For strongly absorbing aerosols, however, s, z) is small and s, z) 

can be as small as 0.09. 


In order to construct model atmospheres sim dating various degrees of contamination 
resulting from polluted air masses or varying weather conditions, we shall do the following: 
(1) use Elterman's values of the total extinction coefficient x(X, z) corresponding to various 
visibility conditions; (2) calculate from the Mie scattering theory the average cross-sections 
for several realistic refractive indices and size distributions, and (3) use these results to 
calculate the single-scattering albedo. 

We can define an aerosol absorptivity parameter f as follows: 

aA\, m, s, z) N(z)c A\, m, s) a A\, m, s) 

“ K^(\, m, s, z) ■ N(z)IT^^^(\, m, s) " ' ' 


Also, since 

tiAa, m, s, z) 

f = 1 “ 

■ K^{\, m, 8, z) 


(44) 


'i 



■1 
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we hiwe 


u)^(X. m. s, z) = 


which reduces to 


z) + (1 - f)#c^(X, m, s, z) 
k(X. V, z) 


<i)^(X. m. s. V. z) = 


£Hj^(X. z) - (1 - z) 


+ 1 - f 


(45) 


(46) 


k(X. V, z) 

Thus, we have the single-scattering albedo for a realistic atmosphere in terms of wavelength 
X, refractive index m, size distribution parameter s, horizontal visual range V, and altitude z. 

It is interesting to consider special cases. If there is no aerosol absorption, = 0 and f = 0 or 

“p(X, z) 

If there is no gaseous absorotion^ -0 and 

K 

z) 


(47) 


u) (X, m, s, V, z) = 1 - f(l - 


V, 2) 


If there is no aerosol scattering, = 0, f = 1 and 

*lp(A, z) 

"o'"'''- " .Toni 

Finally, if there is nogas.;ous scattering, = 0 and 

z) 


w^(X, m, s,V, z) = (1 - f)(l - 


^{X, V, z) 

The Mie scatterir^ computer program was run to obtain values of a 


(48) 


(49) 


(50) 


(X, m, s) and 


_ a, A 

from which the aerosol absorptivity f was calculated. This wat done for 

Deirmendjian’s haze L (continental) and haze M (maritime), as well as for complex refractive 

indices of the following: 


m 




1.5 

1.5 - O.Oli 
1.5 - O.lOi 
1.5 - l.OOi 
1.5 - ooi 


These correspond to all degrees of absorption by aerosols, from no absorption at m = 1.5 to 
complete absorption at m = 1.5 - ooi. For the former, f = 0, and for the latter, f = 1. Figure 7 
Illustrates the spectral dependence of absorptivity for various refractive indices and hazes. 
To get 'exact*' values of cross sections using the Mie scattering program, one would have to 
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run the program intiAitely long. However, sufficiently accurate values can be obtained after a 
reasonable amount of calculation since the values tend to an asymptote. Our program was run 
until the fixed values had only a few percent change as a function of the size parameter x. De- 
tails of the computational procedure can be found in Deirmendjian [10] . 

We can now calculate the single-scattered albedo as a function of wavelength for various 
size distributions, refractive indices, altitudes, and visual ranges. Figure 8 illustrates the 
spectral dependence of for five altitudes. This is for a refractive index of 1.5, that is, no 
aerosol absorption. At the surface, is unity because the ozone contribution is very small. 
However, as we go up in the atmosphere, the ozone absorption band near 0.6 (m becomes 
quite evident. Figure 9 illustrates the same efiect for a heavily contaminated atmosphere with 
high absorption. Here is only about 0.5 near the surface, increasing slightly as we go higher 
into the atmosphere. The low value of near the surface is the result of strong aerosol ab- 
sorption in the dense lower troposphere. 

To see the effects of aerosol and ozone absorption throughout the entire atmosphere, we 
will look at altitude profiles of Figure 10 shows the profile for a hypothetical atmosphere, 
both with no aerosol absorption and with the maximum absorption. Aerosol absorption is 
especially pronounced in the lower troposphere, and the relatively strong ozone absorption 
occurring at higher altitudes is quite evident. 

The profiles for various wavelengths are portrayed in Fig. 11 for a dense haze. This case 
represents the maximum amount of absorption (aerosol and ozone) vhich can occur. It should 
be noted that in the lower troposphere most absorption occurs at the longer wavelengths, since 
aerosol absorption increases with wavelength beyond a certain point (as was seen in Fig. 7), 

At higher altitudes, however, the ozone band near 0,6 jjm dominates. Figure 12 shows the same 
effect except that the case is one of no aerosol absorption. 

We now turn to the more realistic conditions of partial aerosol absorption. Figure 13 il- 
lustrates the altitude profile for three refractive indices: m = 1.5 - O.Oli, or little absorption; 
m = 1.5 - O.lOi, moderate absorption; and m = 1.5 - l.Oi, heavy absorption. Note that there is 
quite a variation in in the lower troposphere. 

The altitude profile for in the first 5 km is shown in Fig. 14 for the two hazes L and M 
and for various amounts of absorption. Varying the size distribution seems to have only a 
minor effect on the profile, as opposed to changes in the imaginary part of the refractive . 

Finally, we can see the effects of a change in visual range on the single scattering albedo 

Figure 15 shows that is essentially constant a& a function of visual range for a given 
haze, increasing rapidly as the amount of aerosol decreases. As any haze (contaminated or 
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FIGURE 13. ALTITUDE PROFILE OF THE SINGLE -SCATTERING ALBEDO 
FOR HAZE L. Visual range = 2 km, complex refractive index m - I 5 - Img. 
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not) decreases, the albedo approaches one char? :terizinjt a pure, aerosol-free Rayleigh 
atmosphere. It is clear that the dominant effect results from a change in refractive index. 

4-5 SCATTERING PHASE FUNCTIONS 

Besides optical depth and the single-scattering albedo, the radiation field within the atmo- 
sphere depends strongly on the single-scattering phase function. Assuming spnerical particles, 
we have used the Mie theory to calculate phase functions for various refractive indices, wave- 
lengths. and size distributions. Figure 16 illustrates the dependence of the phase function on 
the imaginary pari of the refractive index. It is interesting to note that for little or no ab- 
sorption there is the usual peak in the backward direction. For strong absorption, however, 
the backward peak disappears and the distribution becomes more anisotropic. 

The wavelength-dependence of the phase function is exhibited in Fig. 17 for a weakly ab- 
sorbing aerosol. The dependence is as one might expect from elementary scattering theory. 

For shorter wavelengths the function is more peaked than for longer wavelengths. This is 
also true for the backward part of the function. 

The effect of varying particle size distribution can be seen in Fig. 18 for aerosols with no 
absorption. As in the case cf single-scattering albedo, the predominant effect comes not from 
size distribution but from composition. 

Tables of the scattering phase functions for 13 wavelengths, 10 visual ranges, 4 refractive 
indices, 2 hazes, and a complete range of altitudes from 0 to 50 km have been formulated for 
use in computer programs to calculate the radiation field. The implementation of these results 
in our radiative-transfer studies will be presented in Section 5. 

To sum up, we can calculate the important basic optical parameters for an atmosphere 
characterized by scattering and absorbing, spherical aerosol -particles, and by ozone absorption. 
The single-scattering albeiio can be calcul«ited for any refractive index and particle size dis- 
tribution, and we can model any kind of vertically inhomogeneous atmosphere. Likewise, 
single-scattering phase functions can be calculated for any composition and size distribution. 
However, for inhomogeneous, non-spherical particles, actual experimental data on the atten- 
uation coefficients and phase functions must be used. 
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FIGURE 18. SCATTERING PHASE FUNCTIONS AS CALCULATED 
FROM MIE THEORY FOR HAZE L AND HAZE M. Wavelength = 
0.55|i>ni complex refractive index m = 1.5. 
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INHOMOGENEOUS ATMOSPHERES 

We now consider the case of a layered atmosphere in which each layer is homogeneous 
with regard to atmospheric optical parameters. More exact formulations [29) involving per- 
turbation techniques have dealt with the general inhomogeneous Rayleigh atmosphere, but no 
one has solved the general problem for anisotropic scattering. 

5. 1 THE DOUBLE-DELTA APPROXIMATION 

The one-dimensional, radiative -transfer equation can be approximately rolved for hi^ly 
anisotropically scattering atmospheres by using the so-called double-delta approximation 
[30, 31, 32) The resulting equations, however, are different for 1, becoming 

dE’(7) 

= g(T)E;(T) - h(T)E:(T) - Hh{T)E^{r) ( 5 : 


dE’(r) 

- = g(T)E;(T) - h(T)E;(r) + liu^g(r) - H £^( 7 ) (5: 

for the anisotropic field. The functions g(r) and h(7) are 

g(r) = [1- Wq(t)»j(t) ) /Mj, (5‘ 

h(r) = w^(t){ 1 -r/(r)]//i^j (5- 

where rj (r) is an anisotropy parameter. Similarly, the isotropic field E" can be represented 


dE’Mr) 

- 5 ^ = 2pJg(T)E:;(T) - h(r)E"(7)l 
dE”(7) 

=2Mjh(T)E’;(7)-g(r)El'(7)l 


Solving these two sets of coupled differential equations with the appropriate boundary conditions 
allows us tc define the flist "iterative" solution. Knowing the Irradiances, we can form the 
function 
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L(r, /i, [E;(-)6 (m - 

e-(t) + e;;(7) 


(57) 


and insert it into the transfer equation. Representing both the single- scattering albedo 
and the scattering phase function, respectively, by the following step functions: 


’o''’ 

1-1 ^ ^ ^ 


(58^ 


P(7, U. <t>, n', 0') = Vp (/i, 0, /i', 0’)X/_ _ .(r) 

we then can easily solve the radiative -transfer equation for each layer. 


(59) 


5.2 CALCULATIONS 

We now present some results of our exijerience in using the inhomogeneous radiative-ti ar.-^- 
fer equation to find aerosol absorption. Figure 19 shows the spectral dependence of path radi- 
ance and total radiance for two extreme values of refractive indices. Path radiance for the 
contaminated atmosphere is about one-half that for the uncontaminated atmosphere, although 
spectral dependence is nearly the same. Likewise, the total radiance exhibits the same behavior. 

In Fig. 20 we see the spectral dependence of path radiance and total radiance for the same 
conditions as those in Fig. 19, except that the visual range is 2 km instead of 23 km. In this 
case, the effect of aerosol absorption is quite pronounced in that the "clean” atmosphere has 
radiance values many times greater than those for the contaminated atmosphere. 

Figure 21 illustrates the dependence of path radiance on the single-scattering albedo 
Strictly speaking, the curves should coincide at = 1, which implies pure scattering. Because 
of a limiting procedure in the program, however, we used = 0.999. The difference can be 
accounted for m the different phase functions at m = 1.5 and m = 1.5 - l.Oi. Nevertheless, 
the main effect, the strong dependence of path radiance on is quite evident. 

The dependence of both path radiance and total radiance on altitude is exhibited in Fig. 22 
for the two refractive indices For all altitudes both the path radiance and the total radiance 
are greater for the "clean" atmosphere than for the dirty one. 
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FIGURE 21. DEPENDENCE OF PATH RADIANCE ON SINGLE -SCATTERING 
ALBEDO FOR HAZE L. Wavelength = 0.4(im, solar zenith angle = 30®, nadir 
angle = 0°, visual range = 2 km, target reflectance = 0.1, background albedo = 
0.1, complex refractive index m = 1.5 - imn. 
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The same situation is considered in Fig. 23>except that the visual range is 2 km instead of 
23 km. In this case, we note a large difference between the path radiances and total radiances 
for the two refractives indices. Also, one should notice the rapid change in radiance with 
altitude, an easily understandable effect evident in passage through our heavily contaminated 
lower troposphere. Figure 24 shows the same conditions as in Fig. 23, except that we consider 
the variation within the first 5 km. 

Figure 25 illustrates the dependence of path radiance on altitude in the lower troposphere 
for several refractive indices. It is clear that even a small amount of absorption can change 
the path radiance by a significant amount. The results are quJce sensitive to the irs.aginary 
part of the refractive index, especially when 0.01 ^ m^ < 0.10. 

Figure 26 portrays the variations with altitude in path radiance, attenuated radiance, and 
total radiance for the cases of no absorption and heavy absorption. Recall that 

L^ = L^+Lp ( 60 ) 

v;here L^, the attenuated radiance, is the product of surface radiance and transmittance. In 
this case of a heavy haze (that is, for a visual range of 2 km), path radiance increases rapidly 
with altitude for a "clean" atmosphere, and attenuated radiance decreases rapidly. The com- 
bination produces a total radiance which at first decreases, then reaches a minimum, and 
finally increases to its asymptotic limit. In the case of heavy absorption, all radiances are 
lower. 

Finally, Tig, 27 illustrates the behavior of path radiance as a function of the nadir view- 
angle for different refractive indices. The peak occurs at an antisolar angle of 30°. Again, 
contaminated air produces small values of radiance. 

These results clearly show the large differences between radiances for "clean" and "dirty" 
atmospheric conditions. The index of 1.5 represents a silicate haze, not uncommon in various 
parts of the world. The index of 1.5 - l.Oi is quite similar to that for soot particles found in 
the contaminated air of urban areas. Typical refractive indices for various aerosols are pre- 
sented in the work of Kondratyev et al., [33] for the 0.40- to 15.0- jum wavelength region; 
they confirm our original assumption of a very weak spectral variation in the refractive index, 
at least for the visible and near -infrared regions. 

In conclusion, we now have the ability to model any kind of contaminated atmosphere char- 
acterized by aerosol and ozone absorption. The radiation field can vary considerably, depend- 
ing upon the degree of contamination as represented by refractive index and particle size. 
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FIGURE 27. DEPENDEN':E OF PATH RADIANCE ON NADIR VIEW ANGLE 
IN THE SOLAR PLANE R'R SEVERAL REFRACTIVE INDICES. Wave- 
length = O.SSfiin, altitude = 2 km, solar zenith angle = 30°, target reflec- 
tance = 0.1, background reflectance = 0.1, visual range = 2 km. 
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For remote sen.^ing applications, in situ measurements d the air would aid in the definition 
of aerosol content and, hence, would ?»Uow us to calculate the basic radiometric quantities, 
transmittance, irradiance, and radiance. If such mearvurements are not feasible, then we could 
look at multispectral scanner data corresponding to the anti*solar direction in order to obiain 
some knowledge of the phase function. This would help in the estimation of aerosol composi- 
tion. Also, laser scanning or lidar techniques could be employed to determine the same quantity. 
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6 

MULTIDIMENSIONAL RADIATIVE TRANSFER 

In considering the transfer of radiation, we can assume a plane-parallel atmosphere if we 
exclude satellite remote sensing at very high altitudes. Almost all analysis in radiative- 
transfer theory has dealt with a one-dimensional plane-parallel atmosphere, that is, one in 
which the optical parameters, and r, and phase function, p(cos x)> vary only in the direction 
perpendicular to the surface. It has been further assumed that the lower surface Is a perfectly 
diffuse one which is homogeneous in a horizontal plane. These simplifications are sometimes 
necessary in order to solve the radiative -transfer equation. Although assuming that r, and 
phase functions vary only in the vertical direction, we will consider the surface to be a non- 
homogeneous Lambertian type. This set of assumptions will enable us to calculate the influence 
of background reflectances on the target radiance. 


6.1 GENERAL THEORY 

The general, three-dimensional, time-independent, radiative-trajisfer equation for an 
atmosphere wiuch is plane-parallel, scattering and absorbing, but non-emitting is as follows: 

A _ ^ A — _ — A f A A _ A , A 

ft-TL(R. D) = -k(R)L(R, D) + ^ lp(n.fi)L(R, U')d£V (61) 

d’ 

— A 

where Q* is the unit vector in the direction of the scattered photon at position R, and L(R, D) is 
the total spectral radiance, that is, the sum of the direct solar component and the diffuse com- 
ponent. 

The first term on the right-hand side represents a loss of energy along a path, whereas 

A 

the second term represents the sum of all energy scattered into the direction Q. For a one- 
dimensional atmosphere, this can be written as 

L(t. *) - I P(T, Ji. ♦, m’, m’. 

0 -1 



°p(m. 




(62) 


where /i , <f> define the angular coordinates of the scattered radiation 
P(t, fi, (t>, \ 0’) is the single -scattering phase function 

E is the extra-terrestrial solar irradiance of the top of the atmosphere 
^o’ angular coordinates of the sun 

One can immediately find a formal solution to the Integro-differential equation above by 
using an integrating factor; the solution is 
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-(r^-r)/n 

Ljj(7. /j, 0) - Ljj(r^. II. «)e 

2u 1 ~o 

0 -1 

‘V'^o 

+ T — J r e - e e 

4j(p L 

for the upwelliiig radiance. The downwelling radiance is given by 

Lp(r. -p. «) = Ljj(0. -II. 




/i \ 0*)d7* d/i ’ d0* 


(63) 


2 « 1 


J j p(-M, ♦, ii\ p’, 0’)d-'dp'd«' 


0 -1 


‘"o'^O^oP^-'*’ ^ • -Po’ *>0 ^ / -7/w -'K 


Uin - It^) 


- e ]m * li, 


(64) 


There are two general boundary conditions: 
L(0. -/i, 0) = 0 


L(7 


o’ 


2h 1 

/i. 0) - I <^’)[L(T^t 


-/i\ 00 + L (7^, 0O]d/j'd0 

s O 


(65) 


where p*{/i, 0, /i’, 0*) is the bidirectional reflectance of the surface. Thus, Eq. (62), to- 
gether with the bou ndary conditions, can be used to calculate the complete radiation field in 
a plane-parallel atmosphere having horizontal homogeneity. 


6.2 THE UNIFORM DISK PROBLEM 

We can study the effects of background on target by considering a uniform disk with a 
perfectly diffuse reflectance which is surrounded by a Lambertian background surface with a 
different reflectance. First, we shall consider an isotropic scattering law and calculate the 
singly -scattered surface radiance, the singly- scattered solar radiance, the doubly -scattered 
solar radiance, and the attenuated radiance. Secondly, we shall calculate the same quantities 
for an anisotropic scattering law. 
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6.?.1 ISOTROPIC SCATTERING 


Using the one-dimensional solution of Eqs. (63) and (64), we can impose a coordinate sys- 
tem and obtain a solution for the case of a uniform disk surrounded by a spatially infinite, uni- 
form background surface. The appropriate geometry is illustrated in Fig. 28. Let us assume 
that an observer is located at altitude h corresponding to optical depth r and is viewing the 
origin of the coordinate system, i.e., the center of the disk. Radiation emanates from point S 
on the surface, is attenuated as it propagates to point Q, is scattered into direction QP, and 
finally is attenuated to the point of observation at P. Mathematically, we iterate the integral 
equation once for single scattering. Assuming isotropic scattering that corresponds to a very 
light haze condition, we have for the result 


L(7, l)= L(r , l)e 
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(66) 


where 


+ 


E^ and hi are exponential integrals 


The first term represents the radiation from the center of the disk attenuated to point P. The 
second represents the singly -scattered solar radiation. The next two terms represent the 
surface radiation being singly scattered and attenuated to point P. And the last term is the 
doubly -scattered solar radiation. 
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A more advanced case has also been developed for both a finite solid angle of view and 
small scan angles. The correction term to be added to Eq. (66) is too involved to present here, 
but results of the analysis will be. The geometry is illustrated in Fig. 29. The view angle of 
the sensor is Q. steradians, altitude is h (in km), visual range is V, interior Lambertian disk 
reflectance is and exterior Lambertian backgr jund reflectance is In our analysis the 
solid angle U of the sensor is assumea io be 2.5 mrad: the reflectances used are illustrated 
in Fig. 30. !n Fig. 31 we see the singly-scattered surface radiance for black, blue, red, and 
white disks on a green vegetation background. Radiance with the black disk is essentially 
that for green vegetation (as expected); radiance with the blue disk is slightly higher for 
0.4 s X 5 0.65 jum because of the higher reflectance of the disk. For a red panel with a peak 
in the reflectance curve near 0.65/im, the radiance is slightly higher than for a black disk at 
longer wavelengths. For short wavelei^hs (i.e., 0.4 ^ X ^ 0.55 pm) it is essentially zero. 
Finally, for a white disk the radiance is much higher than for a black or colored disk, but it 
generally follows the spectral shape for green vegetation. 

We can now examine the effect of ever-increasing disk size on singly -scattered surface 
radiance. Figure 32 illustrates this effect for a black disk on green vegetation for disk radii 
of 10m, 500m, 1000m, and 5000m. For the smallest disk the radiance is similar to that for 
green vegetation, whereas for a disk of 100 times that radius, the radiance decreases by a fac- 
tor of 8. This Illustrates the strong dependence of surface radiance on disk size. Visual 
range is 2 km. 

Figure 33 is the same situation, except that we now have a white disk on a green vegetation 
background. It should be noted that the effect is just the opposite from that of the black disk. 

As disk size increases, the radiance increases because it is dominated by the high reflectance 
of the white disk relative to the lower reflectance of the background. Also, for small disks the 
curves are similar to those for green vegetation, while for large disks the curves approach 
the spectrally flat white reflectance, biased (of course) by the intrinsic spectral properties of 
the atmosphere. 

These effects are illustrated in a more meaningful way in Fig. 34. Here the singly- 
scattered surface radiance is presented as a function of disk radius for several wavelengths 
by the lower reflectance of the disk. At the peak of red reflectance (X = 0.7 pm), however, 
radiance increases because the disk reflectance is higher than that for green vegetation at that 
wavelength. After a radius of about 5 km, the interference effect hardly changes. 

We can now consider the relative magnitudes of various components of the total spectral 

radiance received by a sensor. Figure 35 illustrates the relative magnitude for a black disk 

on a green vegetation background. L is the singly-scattered surface radiance, L is the 

sur ss 

singly-scattered sky radiance, L^^ is the doubly-scattered sky-radiance, and L^^^ is the su^ 
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FIGURE 28. GEOMETRY FOR ISOTROPIC SCATTERING WITH A 
UNIFORM SURFACE 
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FIGURE 31. SINGLY SCATTERED SURFACE RADIANCES OF COLORED DISKS 
ON GREEN VEGETATION BACKGROUND. Disk radius = 10m, visual range = 2 

km, altitude = 1 km. 
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FIGURE 33. SINGLY SCATTERED SURFACE RADIANCE FOR | 

WHITE DISK WITH GREEN VEGETATION BACKGROUND. Visual . 

rango = 2 km, altltut’ = 1 km. • I 
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of these components. directly attenuated component, is zero in this case because the 

resolution element lies entirely within the black disk. The peak at 3.5/im resuits primarily 
from the peak in extraterrestrial irradiance at the top of the atmosphere. Dominating every- 
thing is the singly-scattered sky radian''e. 

Now we can examine the case of a white disk on a green vegetation background 'see Fig, 36). 
Here is not negligible: in fact, it is the dominant term. The singly-scattered surface radi- 
ance is also not negligible, compared to sky radiance (at least for long wavelengths). 

In all of the analysis so far. only isotropic scattering has been considered — which prob- 
ably tends to overemphasize the surface contribution. In a more advanced treatment we con- 
sider anistropic scattering, as well as various geometries. 


6.2,2 ANISOTROPIC SCATTERING 

We can now consider the more realistic, but also more complicated, case of anisotropic 
scattering. Iterating the formal solution for upwelling radiance produces the following: 
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(67) 


.0 first term on the right-hand side represents the directly attenuated surface radiance* the 
second is the singly scattered solar radiance; the third is the singly- scattered surface radi- 
ance; and the last represents the doubly- scattered solar radiance. We can immediately evalu- 
ate the first and second terms by using actual phase functions; for the last term we can make 
use of the isotropic scattering law. For the third term, however, we must consider an approxi- 
mation for the phase function. 


Let us represent the phase function by the following: 
n-1 -a.//j 

p(/i)=^C^e ' >Cj (68) 
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FIGURE 36. RADIANCE VALUES FOR WHITE DISK WITH GREEN 
VEGETATION BACKGROUND. Disk radius = lOtn, visual range = 

2 km, altitude = 1 km. 
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where C. C^. and a. are constants to be determined. We will attempt to fit the phase function 
in the region 0® to 90^ with the first part of this function and the last half from 90^ lo 180*^ 
with the constant terra C^. This is reasonable since (1) only the first half of the phase function 
is used in the disk problem, and (2) the last half of the phase function is almost constant anyway. 
At the origin, or when ^ = 0^ 

i=l 

Since the phase function is azimuthally symmetric, we have the normalization 


X 

T 2 [ p(p)d/i = 1 


This becomes 




when exponential integral, that is. 




The a. can be determined by using an iteration scheme. The results of this fitting procedure 
are illustrated in Figs. 37-40 for \ - 0.40 and 0.90 pm and for m = 1.5; and m = 1.5 - l.Oi. 

Taking = 0®, or p = 1. we have for the singly-scattered surface radiance: 
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FIGURE 38. EXACT AND APPROXIMATE SCATT 
FUNCTIONS FOR HAZE L AT A WAVE LENGTH 0 
REFRACTIVE INDEX OF 1.5 
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Equation (73) has been used to calculate the singly-scattered surface radiance for combinations 
of spectrally black, gray, and white disks with black, gray, white, and green vegetational back- 
grounds. A black surface is assumed to have a reflectance of 0 at all wavelengths, the gray 
surface a reflectance of 0.5. and the white surface a reflectance of 1.0. The green T^egetation 
background has a reflectance as given iii Fig. 41 [30]. For comparison, Fig. 41 also shows the 
spectral dependence of the extraterrestrial .^olar ir radiance at the top of the atmosphere and 
the product of that ir radiance with the reflectance of green vegetation. The resulting curve 
represents the spectral dependence of the intrinsic radiance for a green vegetation surface 
with no atmosphere. As expected, it is dominated by the spectral character of the vegetation, 
at least for the visible part of the spectrum. 

We will now consider the case of black and white surfaces for varying disk radii. It will 
be assum^'d "^hat the instantaneous field o* view is a circle subtending a half angle of 2.5 mrad 
at the detector. Thus, at an altitude of 1 km the field of view on the surface is a circle with a 
radius of 2.5m. The variation in singly- scattered surface radiance with disk radius is illus- 
trated in Fig. 42 for an atmosphere with a light haze. For the case of a white disk and a black 
background, the sensor ’’sees'' both a small white disk at R = Im and also a larger circular 
ring of black surlace; in addition, the radiance is low. As disk radius increases, radiance 
increases until it reaches a value appropriate to a complete white surface. This occurs at a 
radius of 1 km. The opposite case is also presented, that is when there is a black disk and a 
white background. At small radii the strong influence of the white backgrounc is quite evident, 
even when the disk has a radius of 10m, or four times the field of view. 

The same situation is shown in Fig. 43 for the case of a very ha-^y atmoj-ph^re. Note that 
the radiances here are almost an order of magnitude greater. 

Next, we can see the effect of changing the surface field of view by varying the altitude 
and keeping the disk radius constant at 12.5m. Consider the case of a white disk with a black 
background (see Fig. 44); as the altitude increases the radiance increases, since more surface 
is being observed. However, at an altitude uf o km the sensor "sees” the beginning of a black 
surface and the radiance decreases. At an altitude of 10 km the black area is four times the 
white area and the radiance falls to a very low value. For the case of a black disk and a white 
background, the radiance gradually increases as a function of altitude up to an altitude of 5 km, 
indicating the influence of the white background outside the field of view. Beyond 5 km altitude, 
the sensor "sees ’ a white surface and the radiance Increasei rapidly. 

In Fig. 45, the same situation is represented, except that ihe atmosphere is quite hazy 

An interesting effect becomes evident in this case. As the altitude increases, the radiance 
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FIGURE 42. DEPENDENCE OF SINGI.7 SCATTERED SURFACE RADIANCE ON 
DISK RADIUS. Visual range = 23 km, altitude = 1 km, wavelength = 0.55/im, solar 
zenith angle = 45®, nadir view angle = 0®. 
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FIGURE 43. DEPENDENCE OF SINGLY SCATTERED SURFACE RADIANCE ON 
DISK RADIUS. Visual range = 2 km, altitude = 1 km, wavelength = 0.55fim, solar 
zenith angle = 45°, nadir view angle = 0°. 














disk radius = 12.5m 


2s!i 


willow run laboratories, the university of m cnig^n 


from ;i white disk >icreases up to an altitude of :^bout 0.9 km, after which it gradually decreases. 
This is due to the fact that as altitude increases, more radiation ls included in the acceptance 
cone of the sersor. However, attenuation is also taking place. As a result, after a certain 
altitude more radiation is lost by scattering out of the field of view than is being scattered into 
the field of view. At an altitude of 5 km, for example, the sensor ’’sees” the black surface and 
the rad'anc*.* decreases more rapidly. A similar effect is noted for the case of a black disk with 
a white bac cground. 

We now turn to the cases of black, gray, and white disks with a background of green vege- 
tation. Figun. 46 illustrates the variation ip singly -scattered radiance wdth wavelength for the 
nioe of a black disk. Referring to Fig. 41. we see that the radiance essentially has the spectral 
dependence of the product curv# Co. all disk radii. For radii of 10m and 100m. the *,ctual 
scattering phase function is g. .»er tha the isotropic phase function for all angles appropriate 
to those disks: henc'' the rat'ances are greater. For a disk radius of iOOOm, however, the 
maximum angU ^f sc terin^ at the edge of the disk is 45^. At this angle the actual phase func- 
tion is almost the same as the isotropic phase function; therefore, the radiances for the 'so- 
tropic and the anisotropic cases are almost equal. 

Figure 47 illustrates the same situation, except for a white disk. Here the curves are 
reversed as far as radius is concerned. For large radii the curves have a general behavior 
characteristic of the solar spectrum because the effect of the vegetation background is small; 
for a radius '10m, however, the effect of the vegetatio?. becomes noticeable. The anisotropic 
case * mo pcctral variation than tho isotropic because the actual scattering phase func- 
tiv* ^ ' the f .^rmer are dependent upon wavelength, wlu ?as the isotropic phase function of the 
Ir * .s no spectral dependence. 

Wc can now simulate the case of a sensor in space or at an altitude oi 50 km, which for 
all practical purposes is at the top of the atmosphere. Figure 48 illustrates the spectral ce- 
pender.ee o; -ingly-scattered surface radiance as a function of disk radius for a gray disk. 

Here the effect cf background on target is clear. The field of view i& a circle with a radius of 
^25m. For radii smaller than this, green vegetation as well as the gray disk can be seen; in 
addition, the curve has the spectral dependence of green vegetation. At a radius of 1000m the 
dependence is more closely related to the sc’ar spectrum. 

Figure 49 shows the same effect as for the pi ^vious figure, except that the atmosphere is 
very hazy. In this ca: j the radiance is lower at shorter wavelengths because of the greaUr 
probability of sc ♦♦ering out of the field of view. 

In Fig. 50 we sce the varUition of singly -scattered surface radiance h wavelength for 
black, gray, and white disks, each with a radius of 1000m. For a black disk the only spectral 
varia.tion than can arise is irom the green vegetation background. The rather large 50-perccnt 
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FIGURE 47. DEPENDENCE OF SINGLY SCATTERED SURFACE RADL\NCE 
ON WAVELENGTH FOR A WHITE DISK AND GREEN VEGETATION BACK- 
GROUND. No absorption, solar zenith angle = 45®, nadir view angle = 0®, 
visual range = 2 km, alhtude » 1 km. 
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FIGURE 48. DEPET.'DENCE OF SINGLY SCATTERED SURFACE RADIANCE t 

ON WAVELENGTH AT 50 km ALTITUDE FOR VARIOUS DISK RADT . Visual | 

range = 23 km, no absorption, solar zenith angle - 43^, nadir view angle = 0*^, i 

disk reflectance = 0. 5, bacl^round reflectance = green \ egetation. « 
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FIGURE 49. DEPENDENCE OF SINGLY SCATTERED SURFACE RADIANCE 
ON WAVELENGTH AT 50 km ALTITUDE FOR VARIOUS DISK RADII. Visual 
ra ige = 2 km, no absorption, solar zenith angle - 45°, nadir view angle = 0°, 
disk reflectance = 0.5, background reflectance = green vegetation. 
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gray disk removes this spectral variation and therefore the radiance is essentially the same 
as that for a white disk except for magnitude. 

Figure 51 illustrates the same situation, except for a very hazy atmosphere. Again, the 
radiance at shorter wavelengths is lower due to high attenuation. 

Finally, we can compare the various radiance components (that is, the singly -scattered sky 
radiance, the doubly-scattered sky radiance, the singly-scattered surface radiance, and the 
sum of these three comronents) for a black disk and a green vegetation background (Fig. 52). 
The singly-scattered sky radiance is calculated using the actual scattering phase function; the 
doubly-scattered sky radiance is determined using the isotropic phase function: and the singly- 
scattered surface radiance found via the exponential phase-function approximation. The only 
spectral variation in the sky radiances arises from the solar spectrum. The surface terra is 
small compared to the other radiances, but this is true only because the black disk is so large. 
For smaller disks the surface radiance increases (as was seen in Fig, 46). Thus, the influence 
of outside ref ectance is significant for various disk reflectances and sizes. 

Figure 53 illustrates the same condition, except for a vary hazy atmosphere. In this case, 
the surface term is less important — an effect which is surprising. However, this is for a sen- 
sor at a 50 km-altitude, v;here the sky radiances are larger than at lower altitudes. 

In conclusion, we can say that background surface definitely affects target radiance, the 
amount depending upon the relative size of the target and the field of view. For most cases, 
the anisotropy of scattered radiation is not too important as far as magnitudes of radiance are 
concerned; it is important, however, for the consideration of spectral dependence. 


6.3 THE INFINITE STRIP PROBLEM 

Only in the last few years have investigators tried to find solutions to the multidimensior il 
transport equation. Approximate methods using Fourier transforms have been developed by 
Williams [34], Garrettson and Leonard [35], Krper [36], Er& '"‘nn and Sotoodehnia [37], and 
Rybicki [36). We shall use a similar technique to solve the t -dimensional rrdiative-transfer 
equation for an infinite strip. 


Consider an infinite strip with a reflectance p and a background with reflectance p, as il- 
lustrated in Fig. 54. The position vector R is 

R -- (x, y, z) = R(sin i] cos sin ij sin cos r?) (75) 

A 

and the unit vector is 


il = (SI , ft , n ) = (sin cos 0 , sin 6 sin ^ , cos (^) - (4, p) 
X y z 


( 76 ) 


The radiative-transfer equation becomes 



WAVELENGTH (jum) 


FIGURE 51. DEPENDENCE OF SINGLY SCATTERF D SURFACE RADIANCE 
ON WAVELENGTH FOR VARIOUS DISK REFLECTANCES. Visual range = 'L 
km, no absorption, altitude = 50 km, solar zenith angle = 45°, nadir view angl« 

0°, disk radius = 1000m. 
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FIGURE 52. DEPENDENCE OF RADIANCES ON WAVELENGTH. Visual 
range = 23 km, no absorption, solai* zenith angle = 45®, nadir view angle = 0®, 
disk reflectance = 0, background reflectance = green vegetation, disk radius = 

1000m, altitude = 50 km. 
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WAVELENGTH (|im) 


FIGURE 53. DEPENDENCE OF RADIANCES ON WAVELENGTH. Visual 
range = 2 km, no absorption, solar zenith angle = 45°, nadir view angle = 
0^, disk reflectance = 0, background reflectance = green vegetation, disk 
radius - 1000m, altitude = 50 km. 
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2pl 


“■ T /jLi 

^ z, u, <t>\ t °p{ii, <t>, <l>^) 


2ti 1 


[* p(/i» \ 0')L(x, z, }i\ d<l>' 

6-1 


We can now define a Fourier transfer for the x coordinate. 

X 


L(k, z 
with an Inverse 


ft ♦) ■ f 


L(x, z, Ml 0)dx 


X 

L(x. z, ij, <i>) j e'^r^k, z, M> <^ ) dk 


Taking ths Fourier transform of Eq. (77) leads to the following equations: 

-(t -T-6z)//i 

L(k. z. n, 0) = L(k, 0, /I, i^)e 


.^6(k,ie['-^’^^("-^'’>l/'^L(x,z',M.<Mdz’ 






2ff 1 


"o 0 


(7.) 


(78) 


(79) 


(80) 


1 

I 

'T7 

J 

4 

> 

“I 


11 - 


L(k. z. 4>) = .^6(k) z'. -M. «)dz’ 

z 


2;r 1 X 

I M'.<^')j/l(z')ef^'"’’*'=’^^’’^^|/''L(k, z’. m’. <fr’)dM’d«’dz' (81) 

0-1 z 
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where e = ik| = ik \1 - ^ cos and L(oo, z, #i, 0 ) and L( , z, -/i, <^) are the asymptotic so- 
lutions which coincide with solutions for the one -dimensional problem. By iterating the above 
equations, applying the boundary conditions, and then by taking the inverse Fourier transform, 
one obtains the solution for the two-dimensional equation. This has been done, but the resulting 
mathematical expressions are too involved to be presented here, except for the most interesting 
part dealing with the surface term. Let us consider that term. 

For a Lamoertian surface, L(k, C, u, 0 ) = L(k, 0), and the iterated (first) term beco*::es 
2t: 1 z 

L(k. z. = jp(iu. fi'. «')jf^(z')exp [t - T' !■ c(z - z')]/(i L(k, 0) 

0 0 0 

X ex,. - T' - € ’z' )/(j' d;i' d<(.’ dz' (82) 


The integration over 0 ’ is 
2n 

f exp(€’zV/i’)<i^ ' = 2;rJ 


Fz'kAll - 


(83) 


where is the zero-th order Bessel function. Now, let us represent the scattering phase 


function by 

p(m> 


» -V“' 

e ' 


m'. <!•')= 


(84) 


Kence, 


L(k, z, 1, (*) 

HO 0 




kzjr~2 

M* 


dn’ dz’ 


(85) 


where a = \^ + - t' 

The Integration over p' can easily be accomplished by jslng a Laplace transform; 


, ,, i(k. 0)^^ e*^“ 

L(k, z. 1, = — s— i ) C, \ti(z )e Hu " V ■ 02 

2 


(86) 


It can be shown* that 


*See Appendix A. 
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L(k. O' - 2,7 + |pE(p) - pE(p)]^^ 

Takitig the inverse Fourier transform of Eq. (86) then leads to the following; 



(07) 


z 
C 

where 

M(z', ?) = [k raA|(?/z')^ + lid? (89) 

0 

in which K is the modified Bessel function. The UmHs for M are 
o 

M(z’, 0) = 0; M(z’, co', = ne~'^ (90) 

Therefore, the numerical integration of Eq. (88) will give the spectral radiance for a downward- 
looking observer at any point in the atn^osphere. In particular, we would be able to simulaie 
the response of a sensor as U passes over a natural boundary between tw^o dissimilar surface's 


[M(z’, X + x^) - M(z’, X - dz* (88) 

0 


^ [nE(^) • pE('o)] ^ 
3-1 
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CONCLUSIONS AND RECOMMENDATIONS \ 

i 

As explained in Section 3, Earth's atmosphere has a high degree of variability. In the i 

:^er part of the atmosphere (that is, for altitudes less than 5 km, aerosol density can be very 

high and quite variable in both composition and distribution of particle sizes. Visual range •; 

allows estimation of the aerosol content of the atmosphere but not of composition nor distribu* i 

tion of particle sizes. It is recommended that the atmosphere's spectral optical depth at the alti- f 

! 

bide of the s nsorbe measured and that the spectral q>tical thickness also be measured. This can 

{ 

be done by m ?nitoring solar irradiance at the sensor altitude and the surface. Optical depth is a ; 

. <ore mean* ..^lul parameter than visual range in radiative-transfer studies. In investigations of 

1' ^^e-scale remote sensing by satellites, another quantity which can be used is atmospheric turbidity ! 

.alues of which are becoming available from many measuring sites in the United States. 

! 

In Section 4 we investig...ied the effect of composition and size distribution of aerosols on 
the single-scattering albedo and the single-scattering phase function. Changes in aerosol com- 
position characterized by the imaginary part of the refractive index can considerably alter the 
sin^e-scattering albedo. Of less importance is the change in size distribution. The single- 
scattering albedo is relatively insensitive to visual range up to about 20 km. This is true 
regardless of the degree of contamination (as expressed in terms of the refractive index). 

Similar effects were found for the single-scattering phase functions. The change in size 
distribution did not have much effect on the shape of the phase function, but changes in the 
imaginary part of the refractive index did produce significant differences in the phase function 
at large angles. This suggests that lidar techniques might be employed to measure the back- 
scatter coefficient (or the phase function at 180^). In turn, this would be useful in estimating 
the composition of aerosols. 


In Section 5 we extended the current radiative-transfer model to include absorption. The 
single scattering albedos and phase functions calculated in Section 4 were then used in the 

radiative-transfer equations to determine spectral path radiance and total spectral radiance » 

i 

for various degrees of contamination. Radiances can change by at least a factor of six between > 

a "clean" atmosphere and a strongly absorbing one. Hie actual atmosphere probably lies j 

between these limiting cases. The variation in path radiance with nadir view angle Is not great- | 

i 

ly affected by degree of contamination, but the magnitude can differ by factors of four or five. ) 

h 

i 

In Section 6 we considered the effect of background reflectance on radiance from a target - 1 

element. First, a number of colored disks were chosen as the target with green vegetation as | 

the background. Since this was done assuming isotropic scattering, the spectral variation | 


actually was not realistic. 
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We also considered the more realistic case of anisotropic scattering by fitting actual cal- 
culated phase functions to an approximate exponential formula. The fit was surprisingly good, 
particularly for the angles of importance. These phase functions were then used in the single- 
scattering formulas to calculate singly- scattered surface radiance for a variety of situations. 

We found that bacJcground definitely influences the target and that the amount of influence de- 
pends upon the relative size of the targetand thefield of view, the relative reflectances of target 
and background, and the visual range. Although this influence of surface radiation is usually of 
lesser importance than that caused by radiation from the sky, in the case of large surface 
reflectances the effect is ctf equal importance. The result is that signatures based upon an 
analysis which excluded surface interaction could lead to inaccuracies in the recognition pro- 
cessing of multispectral data. Specific algorithms can be developed, however, to account for 
this effect and correct the data. 

Finally, we investigated the more general problem of determining the radiance at any point 
in tlie atmosphere at which the surface ,s a non-uniform Lambertian reflector. For simplicity 
we assumed an infinite strip with one reflectance and a background with another. Thus, we 
simulated ine case of an aircraft or spacecraft moving across the boundary between two dis- 
similar surfaces. A mathematical expression was obtained for the singly-scattered surface 
radiance as a function a* distance from the strip, altitude of the sensor, reflectances of the 
surfaces, and atmospheric parameters. A numerical integration of the equation is necessary 
for realistic atmospheric conditions. 

We highly recommend that specific flights be made to collect multispectral data for sit- 
uations approximating the geometric and physical conditions used in the analysis of this report. 
Simple experiments will allow us to estimate the degree of contamination, and hence determine 
the variance in radiance. Likewise, the effect of background on target radiance should be inves- 
tigated for the simple geometric conditions we have described. 

Strong gaseous absorption effects in the infrared have been excluded in this analysis, but 
they can be considered by unifying the results obtained here with the ERIM band model for gas- 
eous absorption. The unification of these two models would be particularly advantageous for 
remote sensing applications. Since there is a high degree of variability in the reflective and 
emissive properties of natural materials in the infrared region, a unified model which accounts 
for scattering and absorption by aerosols and gases could be used to correct multispectral data 
over the complete visible and infrared spectral regions. 

In view of the development and use of meterological satellites, we now can obtain meaning- 
ful information on the actual state of the atmosphere at various locations and times over the 
entire globe. If such infoimatlon were readily available to investigator^? they would have a more 
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efficient and economic means of applying atmospheric conditions in the processing of multi- 
spectral data. Therefore, a closer relationship between meteorological investigators and the 
earth-resources remote sensing community is highly desirable. 

As a re.sult of the knowledge and insights gained in this investigation we can now interpret 
real multispectral data in a more meaningful way. By means of techniques developed in this 
report, calculations based upon our present radiative-transfer model can be implemented to 
generate algorithms for the correction erf real data for systematic atmospheric variations. 
Also, variances in multispectral data resulting from atmospheric effects can be simulated in 
a more realistic way. 


102 



FOMMERLY WILLOW RUN LAAORATORtES. THE UNlVCRS'^’^ micmiGAN 



Appendix 

FOURIER TRANSFORM DERIVATION 

The Fourier transform of intrinsic radiance t the surface is derived as follows. 

flC 

L(k, 0)= y^e‘*‘*L(x, 0)dx 

where 0 implies zero altitude. Now, by definition the radiance is 
2 t7 1 

L(x, 0) JjzV(x, 4, 4>f -P’, <^>’)L(x, 0, (f)’)d 4 ’d<^ 

U 0 

which, for a Lambertian surface, becomes 
Hence, 

•Xq X 

_ pE(p)f2 sin kM ^sinkxj] 2 sin kx^ 

= X L k k k — : “ ““ 

= 2x + [pE(p) - 3 e(p)i!^^| 

'i 
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